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■ Abstract 

We study the preparation and manipulation of states involving a small number of interacting 

o\ '■ 

particles. By controlling the splitting and fusing of potential wells, we show how to interconvert 
J> ' Mott-insulator-like and trapped BEC-like states. We also discuss the generation of "Schrodinger 

mi 

cat" states by splitting a microtrap and taking into practical consideration the asymmetry between 

O ; 

the resulting wells. These schemes can be used to perform multiparticle interferometry with neutral 
atoms, where interference effects can be observed only when all the participating particles are 
^ ' measured. 
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Entanglement is at the root of Bell's theorem, which exposes the differences between 
quantum theory and a local classical theory based on elements of reality Q]. The predic- 
tions of quantum mechanics have been experimentally observed with entangled Einstein- 
Podolsky- Rosen (EPR) pairs 0,0] as well as Greenberger-Horne-Zeilinger (GHZ) triples 0. 
A related consequence of entanglement is the possibility of multiparticle interferometry. 
Given a maximally entangled system of iV-particles (a "Schrodinger cat" state) a measure- 
ment of interference between different parts of the wave function corresponding to a single 
particle yields random results. It is only when performing a coincidence measurement on all 



iV particles that an interference pattern is revealed la . Experimental confirmation of this 
result has been obtained using photonic EPR pairs 0, (| and internal states of four ions in 
the same trap but no experiments have been performed using a larger number of parti- 
cles. The latest generation of experiments with photons rely on parametric down-conversion, 
which has the technical disadvantage of an exponentially decreasing number of useful counts 
as N increases. Given that entanglement is the key ingredient in all quantum computation 
and quantum communication schemes, clean experimental studies of its consequences have 
become an active topic of research in the last decade. 

In recent yea, Severn, papers Q Q have suggested the generation of entanglement be- 
tween neutral atoms confined in traps by using their interaction in controlled atomic colli- 
sions. The atoms are guided in their motion and their evolution yields the required entan- 
glement of internal states. Other schemes to achieve this sort of entanglement starting from 

n 

BECs have been suggested [10]. In this letter we present two general iV-atom nonlinear pro- 
cesses. The first one is used to convert a Mott-insulator-like (MI) state into a state with 
all particles in the (many-body) ground state of a single trap (BEC-like state); its reverse 
process converts the BEC state into a MI state. The second process is used to generate a 
Schrodinger cat state starting from a BEC state by controlling the splitting of the well. As 
an application of the processes we discuss a scheme for multiparticle interferometry with 
spatially separated paths. 

In the first process, which is also stage I of the interferometry setup, we start with 
a collection of N atoms in the ground states of N independent traps (MI state). These 
separate atoms can be extracted from a reservoir using a quantum tweezer recently 
proposed by our group. Alternatively, single atoms stored and detected in micro-optical 
traps (which have been experimentally reported but are in excited states of the trap) 
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FIG. 1: Schematics of the multiparticle interferometry procedure. Stage I - creation of N atoms 
in the ground state of the trap starting with N individual atoms in TV" traps. Stage II - creation of 
"Schrodinger cat" state. Stage III - spatial separation of the atoms. Stage IV - applying phases, 
combining on the beamsplitters and measurement. 

can be Raman cooled individually to the ground state. The BEC state (stage I, Fig. [TJ) is 
achieved by bringing together the N wells adiabatically if the interaction between atoms is 
repulsive, as will be shown in detail below. This is a consequence of the quantum adiabatic 
theorem, since the MI state is the ground state when the wells are far apart. Our evolution 
is then represented by 

\w± W2 ■ ■ ■ Wn) — ► = \WW . . . w), (1) 

where the states are properly symmetrized bosonic states. 

In the second process (stage II of multiparticle interferometry) the interaction is switched 
to attractive. This can be done by using a Feshbach resonance pjj. Starting from the BEC 
state, we slowly split the well into two approximately equal mictrotraps, which we label as 
L and R. The lowest energy states are then the ones having all atoms in the left or in the 
right well. Since initially the system is in the ground state, by separating the traps at some 
slow rate v when the wells are far apart we get a linear combination of these two nearly 
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degenerate states, i.e. the system is in the Schrodinger cat state 

-> |^//) =a\LL...L) + Pe i0 \RR...R), (2) 

with a, (3, and 9 real. For perfectly symmetric traps, a = (3 and 6 = 0, but any asymmetry 
makes these parameters rate dependent, as will be discussed in detail below. 

Additional processes are needed to realize multiparticle interferometry. During stage III, 
the interaction is switched back to repulsive and each of the two traps is separated to N. 
This stage can be seen as the inverse of stage I applied to the wells L and R. Again, if the 
separation is done adiabatically the system remains in the ground state which in this case 
corresponds to a single atom in each one of the wells. The state is now 

|*//> -> \*m) = a \L X L % ...L N )+ f3e ld \R X R 2 ...R N ). (3) 

Subsequently, the atoms in wells derived from the original R well are subjected to additional 
phase shifts <j>i, 2 • • • <Pn, which can be applied, for example, by adjusting the depth of the 
wells adiabatically. 

In the final, stage IV of the scheme, we combine states and Ri in a 50-50 beamsplit- 
terQ. Notice tteint fce e^e* 0n l y one of these two is occupied so the interatomic 
interaction plays no role in this stage. We denote the outputs of each beamsplitter by A{ 
and Bi and assign a value of +1 to the measurement of an atom in channel A iy and —1 to the 
measurement of atom in channel Bi. The probability, P(+l), that the product of all measure- 
ments gives +1 (for instance A1B2B3 in the case of three atoms) is (1 — a/3 cos (A + 6)) /2, 

TV 

where A = 0j. The probability for the product to be —1 is P(— 1) = 1 — P(+l), hence the 

i=l 

expectation value over a large number of measurements is —a/3 cos (A + 9). We would like 
to stress that a correlated measurement of less that N atoms does not show any dependence 
on phase and appears random. 

In order to obtain the relevant parameters for the operation of our inteferometer, we 
study the evolution of an N particle system using optical microtraps. As an example, we 
numerically solve the Schrodinger equation in the case of three atoms in a quasi-lD config- 
uration. This is achieved by strongly trapping the atoms in the perpendicular dimensions, 
effectively freezing these degrees of freedom. We scale the equations choosing units of length 
L u = 2 /Jin, of energy E u = h 2 / (2M u Ll) and of time t u = h/E u . The particle interaction is 
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FIG. 2: Stage I: adiabatic energy levels for three atoms in three wells with repulsive interaction 
as a function of d. The other parameters of the potential are V$ = 10, a = 0.5, Uq = 10, 
<73 = -Qi = 10~ 4 , <72 = 0. 

represented by a delta-function potential 

U(x 1 ,x 2 ) = U 5(x 1 -x 2 ). (4) 

The atoms are also subject to external potentials due to the optical traps, which in each 
stage are 

V Ijn (x, d) = £ (1 + qi )V(x, (i - 2)d), 

i=i (5) 

V n (x, d) = (l + qi )V(x, -d/2) + (1 + q 2 )V(x, +d/2), 

with 

V(x,d) = -V ev^-^±^j. (6) 

The qi parametrize the asymmetry between the intensities of the beams denning the different 
wells; we assume that these are 10~ 4 . 

Let us consider first the evolution during the first and the third stages of the operation. 
There are four different energy scales in the problem. The first one is the energy difference 
between the energy levels localized in different wells, which we can estimate as E asym qV . 
The second one is the energy required to move one of the atoms to an already occupied well, 
estimated to be E- mt U /a where a = (Vo/a 2 ) 1 ^ is the width of the wave function in a 
well. The third scale is the energy E^c ~ a o 2 required to put one of the atoms in an excited 
state of one of the traps. The last energy scale (E D ~ (ir/ND) 2 ) is the energy required to 



excite the atoms out of the ground state when the distance between the wells is D »s 2a, at 
which time the trap can be approximated by a square well of width ND. We operate in the 
regime in which 

Ezsym "C -Ejnt; -^exc; Ej). (7) 

Figure El shows the dependence of the adiabatic levels on the separation d during this stage. 
The presence of the small asymmetry in this stage does not affect the nature of the ground 
state, which is non- degenerate. Joining or separating the wells at a slow speed keeps the sys- 
tem in the ground state, i.e. the lowest curve in the figure. We can estimate the rate at which 

n 

the adiabaticity is lost by applying the Landau-Zener formula [16], v a d ~ (AE glL p) 2 / (dE/dx). 
The slope can be estimated as y/NVo/cx 2 . The size of the gap depends on which of the three 
large energy scales in (JZj) is the smallest. In the example that we are presenting, all three are 
roughly the same order of magnitude. The probabilities |a;| 2 = | (^l^) | 2 to find the system 
in the states \ipi) at the end of the evolution is plotted in Figure 01 as a function of the speed 
v. In our example, the critical rate is v c i^ = 0.35; the probability to find the system in other 
states is less than 0.01. For multiparticle interferometry it is critical not to accumulate an 
additional phase during the third process due to the asymmetry between the right and left 
set of wells. This gives rise to a lower bound for the allowed velocity, as explained below. 
For the parameters chosen in the figure this is v c ^i = 0.09. 

Between these stages and stage I we need to change the sign of the effective interaction 
between the particles. For the cases we are considering, the particles remain the ground state 
with very high probability (of the order of 99%) even if this change is performed suddenly. 

During stage II the adiabatic energy levels as a function of d are shown in Figure 0] 
Once again we have four energy scales, which can be approximated by -E a sym ~ NqVo, 
E- m t ~ (N — l)\Uo\/(Jo, E exc « <Jq 2 , and Ed « (tt/2D) 2 . Once again, we work in the regime 
in which (J2J) is valid. Separating the wells adiabatically maintains the system in the ground 
state, which corresponds to all N atoms being in the lowest of the two wells, which is not 
the desired state. In order to mix the lowest two energy states we need to evolve the system 
non-adiabatically with respect to the lowest gap but at a slow enough speed to remain 
adiabatic with respect to the larger gap. Below v c n,2 = 0.27 the probability to tunnel to 
these excited states is less then 0.01 and entanglement is obtained with a/3 = 0.99 or larger. 
On the other hand, the asymmetry yields a dephasing between the two parts of the wave 
function 6 = -E asym t S ep, where the separation time is inversely proportional to the velocity v. 
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FIG. 3: Stage III: probabilities to find the system in the adiabatic states after a single well with 
three atoms is split into three wells with an atom per well (dfmai = 3.0) as a function of the speed 
v. The energy levels are the ones shown in Fig. [3 For velocities smaller than denoted with dashed 
line probability to state in the ground state is larger then 0.99, for velocities larger than denoted 
with dashed-dotted line dephasing is less then 0.1. For stage I the dynamics are very similar except 
there is no limit on how slow the process could be done. 

Allowing a maximum dephasing max , we must go faster than v c i t \ ~ qVo } niN/(f) milx . This 
calculation assumes, however, that the asymmetry is constant. In a practical situation, q is 
driven by fluctuations in the laser power, and consequently the phase 9 grows diffusively, as 
the square root of t sep instead of linearly, making the condition less restrictive. 

The only two conditions for the applicability of the method are related to the asymmetry 
of the potential. As long as condition (|7J) is met and as long as v c ^ is larger than t> Cjl , 
there is a range of velocities for which the operation is successful. The critical velocities 
have different dependence on N, so for fixed values of the parameters defining the potential 
and the interaction, there is a largest number of atoms for which this happens. However, 
by choosing a different set of parameters this condition can be relaxed. In particular, the 
strong N~ 2 dependence of the preparation of the MI state can be overcome by separating 
the atoms in series instead of doing it parallel (for N = 2 n , we can think of n steps in which 
each well is split into two). 

Finally, numerical values for realistic experimental parameters are given. In the model 
described above the effective interaction between atoms is determined by the scattering 
length a and the strength of the confinement in the transverse direction. The frequency u± of 
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the confinement in the case in which the system stays in the ground state of transverse motion 
may be expressed in terms of the dimensionless interaction parameter Uq used above as 

U± 4\a\M u L u 1 ' 

Hence it is desirable to use atoms with the largest product of mass and scattering length 
possible. In Table |J we present the rescaled values used in the calculation for two workhorses 
of cold atom experiments, sodium and rubidium. The magnetic fields needed to observe 
Feshbach resonances in alkali atoms are typically hundreds of gauss Q|. In the proposed 
scheme for multiparticle interferometry one should work on the side of the resonance where 
the scattering length changes sign to avoid the losses associated with crossing the resonance. 

We acknowledge support from the NSF, the R.A. Welch Foundation, MGR also acknowl- 
edges support from the Sid W. Richardson Foundation. 
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FIG. 4: Stage II: adiabatic levels of three atoms in two wells in the case of attractive interaction 
for different values of the separation d . The other parameters of the potential are Vo = 30, a = 0.5, 
U Q = -4, Ql = 0,q 2 = 10- 4 . 
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FIG. 5: Stage II: full lines are probabilities to find the system in the adiabatic states after the 
separation of one well with three atoms to two (dfi na i = 3.0) as a function of the speed v. The 
dashed line is 9. The interaction is attractive and the parameters are the ones used in Fig. For 
velocities in the interval between vertical lines the desired state is prepared with probability of 0.99 
and dephasing smaller than 0.1. 
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